Intrinsic self-rotation of BEC Bloch wave packets 
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The semiclassical theory of Bloch wave packet dynamics predicts a self-rotation angular mo- 
mentum in asymmetric periodic potentials, which has never been observed. We show how this is 
manifested in Bose-Einstein condensed atoms in optical lattices. Displacing the wave packet to a 
corner of the Brillouin zone we obtain a current distribution with a non-quantized angular momen- 
tum, independent of the size of the distribution. A weak interatomic interaction does not modify 
the results, affecting only the rate of spreading in the lattice. A strong repulsive interaction results 
in a collapse of the wavefunction into matter-wave lattice solitons. 



The semiclassical equations of motion for electrons in 
crystals have played a fundamental role in the theory of 
charge transport in metals and semiconductors Be- 
cause of Berry phase and gradient corrections UGl , two 
striking effects occur for systems without time reversal 
or spatial inversion symmetry. The first one is a lat- 
eral displacement of the particles relative to the external 
force, and is responsible for the anomalous Hall effect in 
ferromagnets |j, |5|. The second one is the appearance 
of an orbital magnetization, due to the rotation of the 
wave packet around its center. This latter effect has re- 
ceived little attention and no experimental observation 
of it has been performed so far. In this letter we propose 
a demonstration of this wave packet self-rotation using 
ultracold atoms in optical lattices. Such a system has 
been employed in the past to exhibit a number of basic 
phenomena in condensed matter physics |g, LZj • 

A dual purpose of this work is to study the generation 
of novel matter-wave states using the knowledge gath- 
ered in electron transport theory. Bose-Einstein conden- 
sates can be manipulated experimentally with great pre- 
cision and the equation of motion for the condensate 
wavefunction in a periodic potential coincides with the 
one for electrons moving in a crystal. Several experimen- 
tal studies of condensates in optical lattices have already 
been performed |3, llOj . The generation of states in an 
asymmetrical potential must take into consideration the 
two effects mentioned before. 

The semiclassical description assumes that the wave 
function for a particle is a superposition of Bloch waves 
from a single band 
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in which the distribution /(k) is centered at a point k c 
in the Brillouin zone with a small dispersion tTfc. This in 
turn yields a spatial dispersion a r much larger than the 
size of a unit cell. The Bloch waves have the property 
that | = e lk ' r |uk), where Uk(r) preserves the periodic- 
ity of the lattice potential Vj a tt (r) . These wave functions 
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We restrict our study to a single band, thus dropping the 
band index n in what follows. 

By defining the Lagrangian £(k c , r c , k c , r c ) = 
(^f\(ih-^ — H)\^) using as coordinates the position of the 
center of the wave packet in real and reciprocal space, we 
can find the equations of motion, which read 
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For atoms in an optical lattice, the "electric" and "mag- 
netic" forces are the inertial Coriolis forces, if we describe 
the motion in a frame moving with the lattice which is lin- 
early accelerated and/or rotated, respectively. Equation 
Q is the Lorenz force equation found in any solid-state 
physics textbook, but J2J) contains two corrections from 
the standard result. One is the presence of the Berry 
curvature, defined as 
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(notice that this can be nonzero for complex u(r)). It 
acts dynamically as the symmetrical analogue of the mag- 
netic field if we exchange the r c and k c variables, and is 
related to the Hall effect in magnetic sub-bands Q and 
the anomalous Hall effect in ferromagnets 0, 0|- The 
second correction is to the band energy, which contains 
a term proportional to the orbital magnetization energy, 
£s(k c ) — £(k c ) + 2lfS(k c ) • B(r c ). This magnetization 
is proportional to the internal angular momentum of the 
wave packet, i.e. 



S = 



J dr (r - r c ) x J(r), 
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where J(r) is the current density 0- In terms of the 
Bloch waves, it is given by 

s(kc) = *y { wJ x {H{kc) £{kc))l W c } - (7) 

As long as the semiclassical approximation holds, this 
result is independent of the distribution /(k) used, and 
hence of the width of the wave packet in real space. In 
particular it is independent of time for a fixed value of 
k c ; the wave packet maintains its angular momentum 
constant as it spreads on the lattice. 

The reason that in most circumstances these correc- 
tions to the semiclassical equations of motion do not need 
to be taken into consideration is the constraints imposed 
on them by the symmetry of the Hamiltonian. If the 
system possesses time-reversal symmetry, then both vec- 
tors f2 and S are odd functions of k c . On the other 
hand, if it possesses inversion symmetry then they are 
even functions of k c . This implies that in systems with 
both symmetries they vanish throughout the Brillouin 
zone. 

We will consider a system of Bose-condensed atoms 
moving in an asymmetric two-dimensional optical lattice. 
The interaction between the atoms is neglected for the 
moment. In order to break inversion symmetry, we con- 
sider an asymmetric hexagonal lattice, as shown in Fig. 
^a). In the tight-binding regime, in which the poten- 
tial wells denoted by blue color are deep, the atoms are 
located in the ground state of these wells. Since the cal- 
culations can be performed analytically in this case, we 
consider it first jllj . 

Let the on-site energies for sites A and B be ea — e g /2 
and 6b = — £g/2, and the tunnelling matrix element be 
non- vanishing only between nearest neighbors with mag- 
nitude h. The lattice is triangular, with basis vectors 
ai = \Z3ae y and a.2 — |ae x + ^-ae y ; here a is the dis- 
tance between nearest neighbors. Numbering the sites in 
the lattice by (to, n) such that R TOi „ = mai + n&2, we 
obtain the discrete Schrodinger equation for the eigen- 
vectors of the Hamiltonian 

£ i>i,n = *A 4>i, n ~ ft(Vm,n+l + + <+l,J> ( 8 ) 
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The Bloch waves for this Hamiltonian are calculated 
by writing ipmn = <i> A ' B exp(ik • R m .„), with which we 
obtain 
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The off-diagonal element is V(k) = -/i(l+e lk ai +e lk a2 ). 
The diagonalization of the Bloch Hamiltonian is straight- 
forward, yielding the two energy bands 
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FIG. 1: (a) potential energy seen by the atoms; the blue 
color corresponds to the lowest values. In the tight-binding 
approximation, particles are localized in the A and B sites 
only. The corresponding Brillouin zone in reciprocal space is 
shown in (b); due to symmetry, all points marked by a circle 
correspond to the same K point. In (c) we show the band 
structure (solid line, left axis) and the angular momentum S 
(dashed line, right axis) for a tight binding band with e g = 
h = 1. 



The energy bands are depicted in Fig. ^c). The bot- 
tom band has a minimum at the origin (r point) and a 
maximum at the K point, ko = 47r/(3-\/3a) e y , at which 
T^(ko) = 0. This is also the point of closest approach to 
the excited band, with the gap being equal to e g . 
For the ground state wave function, 
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In two dimensional systems, the Berry curvature and the 
angular momentum are directed along the perpendicular, 
z-direction. Their values for the system are 
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5(k) = — (2f+)n(k). 

An intriguing relation at the K point, where S is mini- 
mum, is 
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where M* is the effective mass, given by M* = 
-2e g h 2 /{9h 2 a 2 ). 

The angular momentum S is plotted in reciprocal space 
in figure ^c) (dashed green line). From (|14f> we see that 
the dependence of the Berry curvature on quasimomen- 
tum in reciprocal space is very similar. In particular, 
both of these quantities attain their maximum absolute 
value at ko and all symmetrically located points. The 
value attained diverges as the gap size e g goes to zero. 
Notice that the angular momentum is not quantized. 

The angular momentum that we have calculated is the 
total value for the wave packet. To study its distribu- 
tion in real space, we consider a Gaussian distribution of 
Bloch waves around the point ko, i.e. 
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In figure |3 a) we show the distribution of currents as a 
function of position, which show a clear rotating pat- 
tern. This distribution of currents is not the one found 
near a vortex, for which the currents diverge at the cen- 
ter. Approximating the off-diagonal term in the Bloch 
Hamiltonian near ko as 

V(k) a i^V*/6 ((fc x _ ko , x ) + i(k y - k , y )), (17) 

we can calculate the wave function coefficients as 
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Notice that the probability for B sites has a maximum 
at the center of the wave packet, while the probability 
for A sites is zero there. This is related to the fact that 
for a single Bloch wave at ko the A sites are completely 
depopulated, and it is only through the presence of a dis- 
tribution of Bloch waves that some of the A sites acquire 
some density. The distribution of A sites has a maxi- 
mum when R m ,n ~ l/cfc- Since the currents depend on 
the values at neighboring sites, it is at this distance that 
the currents are maximized. 

In a real experiment, the atoms will move in an optical 
lattice generated by counterpropagating laser beams or 
a holographic method, in which the potential is gener- 
ated with a phase or amplitude mask 12]. The potential 
to be chosen can be obtained by calculating the Fourier 
transform of a sum of Gaussian wells located at the A 
and B sites and keeping the largest terms. The potential 
obtained is of the form 

Matt(r) = Vb[cos((k a - k 2 ) • r) + cos((2k 1 + k 2 ) • r) 
+ cos((2k 2 +ki) ■ r)] 
+Vi[cos(ki • r + a) + cos(k 2 ■ r + a) 
+ cos((k 1 +k 2 ) -r-a)]. (20) 




FIG. 2: (a): distribution of currents in the tight-binding limit 
for a wave packet centered at ko. (b): calculation of the 
averaged angular momentum distribution for atoms in the 
periodic potential H2UH with Vb = —4, Vi = 4.664 after the 
wave packet has been displaced to ko. The averaging assumes 
an experimental precision in position of CTinst = a. 



Here the wavevectors are ki = — l^e-r + -?2-e,,, and k 2 = 

|^-e x . We will use a system of units in which % = M = 
a = 1 . It is worth noticing that the ratio between the two 
amplitudes V\ and Vb controls the widths of the gaussian 
potentials located at the lattice sites, while the phase 
yields the ratio between the depths of the potentials. For 
the symmetric potential, a — In our calculations 

we have used a value of a = 1.1. 

A wave packet of the form Q at k c = in the low- 
est energy band of a periodic potential can be obtained 
starting with a Gaussian wave packet in real space (e.g. 
the ground state of a harmonic trap) and adiabatically 
turning on the lattice . The study of the angular mo- 
mentum as a function of the quasimomentum of the wave 
packet can be performed by applying an "electric" field 
(see Eq. (0J), which in the case of cold atoms corresponds 
to an acceleration of the lattice. The atoms are displaced 
in reciprocal space along the direction of the acceleration, 
while they perform Bloch oscillations in real space. 

The calculated motion of the wave packet in real space 
and its angular momentum as a function of time are 
shown in Fig. [3] The wave packet starts at I, moving 
against the direction of the acceleration. As the quasi- 
momentum reaches the K point in reciprocal space (II in 
real space), the Berry curvature displaces the distribu- 
tion to the right, and the angular momentum achieves its 
smallest value. The averaged distribution of the angular 
momentum is shown in the right panel of Fig. [21 which 
agrees with the result for the tight-binding model (left 
panel of the same figure). Eventually, the wave packet 
reaches point III, where it shows no rotation and starts 
retracing its path. When it reaches II again, it possesses a 
maximum (positive) angular momentum, which it looses 
once it comes back to I. 

The velocity distribution of the atoms (from which 
the angular momentum of the wave packet can be cal- 
culated) can be measured by Doppler-sensitive Raman 
transitions In this way, atoms with a selected veloc- 
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FIG. 3: (a): motion in real space of the center of the wave 
packet, (b): angular momentum as a function of the position 
in reciprocal space, calculated for an "electric" field of — eE = 
0.05 applied along the y-direction. 

ity are driven to a dark state, and one can get rid of the 
rest of the atoms by shining resonant light. Imaging the 
remaining atoms, one may obtain the real space distri- 
bution of atoms with the given velocity. By changing the 
detuning, different sets of atoms are observed. Another 
observable result in this case is the motion caused by the 
anomalous velocity (due to the nonvanishing Berry cur- 
vature). This displaces the center of the wave packet by 
an amount equal to 

Ar c = - [ dk' x n(k'), (21) 
Jo 

an amount independent of the magnitude of the driving 
-eE field. 

Care must be taken in the experiment not to excite 
the atoms to a higher energy band while the lattice is 
accelerated. The largest probability for this to happen 
occurs when the wave packet is at ko, where the gap 
between the two bands is smallest and the particles can 
undergo Landau-Zener tunnelling. The probability for 
this to happen is proportional to exp(— e 2 /\eE\), there- 
fore the electric field must satisfy \eE\ <C e 2 ,. 

One has to be careful in the design of the lattice in 
which the rotation is observed in order to make a di- 
rect comparison with the tight-binding model. Using a 
potential depth Vq much smaller than the one we used 
in the simulations one has to include the probability of 
atoms tunnelling to a non-nearest neighbor. These cur- 
rents between B sites would not add to the total angular 
momentum in the wave packet but can yield a much more 
complicated distribution of currents. 

We finally comment on the consequences of adding the 
interatomic interaction between the particles. The inter- 
action corresponds to the addition of terms g\ r ^ m ,n\ 2 '4'm,n 
in the right hand sides of © and ©■ For weak inter- 
action strength we observe no change in the values of 
the angular momentum at kg , which remains constant in 
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FIG. 4: (a) Dispersion of the wave packet with time for dif- 
ferent values of the interaction strength in the tight-binding 
limit (e s = 6, t = 1). The wave packet is accelerated from 
the r point to the K point, where the interaction strength is 
adiabatically increased. The lines correspond to g = (blue 
line), g = — 1 (red line), g = 1 (green line), g — 1.9 (cyan line), 
and g — 3 (magenta line). In this last case, part of the wave 
packet collapses to a soliton when the size of the wave packet 
is minimum, with the remaining part diffusing at a fast rate 
afterwards, (b)-(e) show the wave function in real space after 
dispersing in the lattice for t « 200; the interaction strengths 
are g = 1.5 (b), g — 7.5 (c), g = 15 (d), and g = 30 (e). 



time. The rate of dispersion of the wave packet is mod- 
ified, with a repulsive interaction (<? > 0) slowing down 
the dispersion and an attractive interaction (g < 0) ac- 
celerating it (see Fig. 0J. This counterintuitive result 
is related to the negative effective mass associated with 
particles near the top of the energy band 0|- Moreover, 
for repulsive interactions larger than a critical value, the 
system exhibits collapse into one or more breathing soli- 
tons ^(|. These can be related to solitons in the gap 
between the two bands ^3] similar to recently observed 
two dimensional optical solitons 0] and will be consid- 
ered in a separate article. 
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